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ABSTRACT 

It has been proposed that the excess in cosmic ray positrons at approximately 8 GeV, 
observed on both flights of the HEAT balloon experiment, may be associated with the 
annihilation of dark matter within the Milky Way halo. In this paper we demonstrate 
how the self-annihilation of neutralino dark matter within local substructure can ac- 
count for this excess, and estimate the annihilation cross-section for several benchmark 
minimal supcrsymmctric (MSSM) models. We also demonstrate the effect on the per- 
mitted parameter space as tidal stripping effects and destruction of substructure by 
mergers becomes increasingly severe. 
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1 INTRODUCTION 

A large body of evidence pertaining to the existence of non- 
baryonic cold dark matter (CDM) has been established over 
the past several de cades, including the la rge-scale distri- 
bution of galaxies dAbazaiian et all l|2003l)), the study of 



primordial light element abundances jCvb urtet ahl (12004 
galact ic rotation curves dBe geman et al supernova 



data dPerlmutter et al. (1999)) and the power spectrum of 
anisotro pies observed in the cosmic microwave b ackground 
(CMB) dBennett et al.l d2003Mde Bernardis et ail d2000|) ). 
At the 2a confidence level, the CD M energy density para m- 
eter is now ^cdm/i 2 = O.liat^ois dBennett et alJ d2003ft ). 

A popular candidate for CDM is the lightest super- 
symmetric (SUSY) neutralino, which is a superposition of 
Higgsinos, Winos and Binos. Consequently the neutralino is 
electrically neutral and colourless, only interacting weakly 
and gravitationally and therefore very difficult to detect di- 
rectly. In R-parity conserving supersymmetric models the 
lightest n eutralino, being t he ligh test SUSY particle (LSP), 
is stable dWeinberg et alJ d2000l) ). Consequently, in a sce- 
nario where present-day CDM exists as a result of thermal- 
freeze out, the dominant species of CDM is likely to include 
the LSP. The relic density of the LSP will then heavily de- 
pend on its mass and annihilation cross-section. The neu- 
tralino is a popular candidate for CDM because the expected 
values of these parameters are such that the corresponding 
relic densities are of the same order as the currently accepted 
value of ficDM- 

In 1994 and 1995, the High-Energy Antimatter Tele- 
scope (HEAT) observed a flux in cosmic ray positrons, 
very much in excess of theoretica l predi ctions, and peak- 
mg 



close to 8 GeV dCoutu et al] dl999ft ). This observa- 



tion was con fi rmed by a second flight by HEAT in 2000 
dCoutu et al.l d200l]) ' ). Th e effect was originally predicted 
in 1984 TsluTetHnT (I1984I) ). for both cosmic ray positrons 
and antiprotons. While the cosmic ray antiproton feature 
remains elusive, at least for theorists it is now a well estab- 
lished idea that local dark matter annihilations are a pos- 
sible source of the positrons contributing to the observed 
excess. Of course, improved data is urgently required to val- 
idate the significance of this observation, and these are ex- 
pected to be taken in the near future by experiments like 
PAMELA, AMS-02 and Bess Polar. In this paper, we sum- 
marise why extant models for local dark matter annihila- 
tions fail to adequately account for the current spectral data 
on the positron feature, and present a new model. 



Studies in which the local (within several kpc) distri- 
bution of dark matter is considered to be smooth have con- 
cluded that the dark matter annihilation rate is insufficient 
to reproduce the observed positron flux. This naturally leads 
to the idea of a non-uniform, or "clumpy" , dark matter dis- 
tribution providing the overdensities necessary to elevate 
dark matter annihil ation rates su fficiently to produce the ob- 
served excess dBaltz et al] d200 1]) ). The distribution of such 
substructure can then be estimated through the use of simu- 
lations. Unfortunately, while simulations generally agree on 
the distribution of dark matter clumps, there is a divergence 
of opinion on the survival of the lightest clumps after the ef- 
fects of hierach ial structure formation and tidal strippin g 
are considered i Diemand et al] d2005T) . IZhao et all |2005), 
iDiemand etafl fcOOd) ). 



* E-mail: dtc@astro.ox.ac.uk 
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We note that throughout this paper we shall consider 
the LSP to be the lightst neutralino and therefore use the 
two terms synonymously. 
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Figure 1. The positron energy spectra from neutralino annihi- 
lations for modes utilised in this paper. Solid lines represent the 
positron spectra, per annihilation, for \\ —* bb, for neutralinos 
with masses of 50 , 150 and 600 GeV. Dotted lines are the same, 
but correspond to the process \X ~ * t + t - . Dashed lines cor- 
repond to the positron spectra fom the process \X ~ * W~^~W~ , 
for a 150 GeV neutralino, where the spectrum from the process 
XX ~~ * ZZ is very similar. 



2 POSITRON SPECTRA FROM 
NEUTRALINO ANNIHILATION 

Positrons can be produced in several neutralino annihila- 
tion modes. For example, monoenergetic positrons can result 
from the decay of gauge bosons produced in the processes 
XX — * ZZ or XX ~ * W + W~ , producing positrons of energy 
m x /2, where m x is the neutralino mass. A small amount of 
positrons can result from the direct channel, XX e + e~, 
however such a channel seldom makes a significant impact 
on the overall positron flux. A continuum of positrons, ex- 
tending to much lower energies, can also be produced in the 
cascades of particles produced in annihilations. The spec- 
trum of positrons produced in neutralino annihilations can 
vary significantly depending on the mass and annihilation 
modes of the LSP. 

If the neutralino is lighter than the W ± and Z bosons, 
annihilations will be dominated by the process XX ~ * bb 
with a minor contribution by XX ~* t + t _ . Assuming anni- 
hilations are dominated by the former process, the result- 
ing positron spectrum will depend entirely on the mass of 
the LSP. For heavier LSPs, the annihilation products be- 
come more complex, often determined by several dominant 
annihilation modes including XX ~> W + W~ , XX ~* ZZ 
or XX ~ * as weu as XX — * bb and XX ~ * t + t~ . In 
our calculations, the positron spectrum per annihilation 
per energy inter v al d<b /dE R +, uses results from PYTHIA 
ilSiostrand et alJ fcOOlIf ) . imp lemented in the DarkSUSY 
package iG^ndoloetaLl fcOOfj) '). 

In figure^ we display results for positron spectra from 
neutralino annihilatio ns for the most impor tant annihilation 
modes as produced bv lHooper et alJ (I2004T) . Solid lines rep- 
resent the positron spectrum, per annihilation, for XX ~ * bb, 
for LSPs with masses of 50, 150 and 600 GeV. The dot- 
ted lines are the same, but from the process XX ~ * t + t~. 
Gaugino-like annihilations typically produce a spectrum 
which is dominated by bb at low energies, with contribu- 
tions from t + t~ only becoming important at energies above 
about half the LSP mass. For neutralinos with a dominant 
higgsino component, annihilations to gauge bosons often 



dominate. Dashed lines represent positron spectra resulting 
from the process XX W + W~ for a 150 GeV neutralino, 
where the spectrum corresponding to the process XX ~* ZZ 
is very similar. 



3 POSITRON PROPAGATION MODEL 

Cosmic ray positrons being electrically charged particles 
diffuse through the interstellar medium (ISM) under the 
electromagnetic influence of the galactic magnetic field and 
background radiation. The spatial distribution of the galac- 
tic magnetic field is extremely complex, resulting in positron 
trajectories which are well approximated to random walks. 
Whilst traversing the ISM, the positrons radiate through 
synchrotron losses owing, to the surrounding magnetic fields, 
and by inverse Compt on scattering off of sta rlight and CMB 
background radiation JWebber et al] <1992h L 

The diffusion-loss equation describing the above process 
is given by 



d dn e + 
didE„+ 
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where dn e + /dE e + is the number density of positrons per 
unit energy interval, K(E e +,r) is a diffusion coefficient, 
b(E e +,r) is the rate of energy loss and Q(E e +,r) is the 
source term. 

We parameterize the diffusion coefficient and the rate 
of energy loss as follows 
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based upon measurements of stable nuclei in cosmic rays 
(prima rily by fitting to obs ervations of the boron to carbon 
ratio) ijMaurin et alJ l)2002h 'l. where in equations (2) an d © 
we have assumed a spatially uniform galactic magnetic field. 
We then define K and b to be constant within a "diffusion 
zone" , which here is considered to be a radially-infinite cylin- 
drical slab of thickness 2L ~4kp c, which is the best-fit to 
observations iMaurin et alJ fc002lLlWebbe 



aer et alJ J1992 
We solve equation Q for the local positron flux, using 
the proceduce described in the appendix to this paper. 



4 POSITRON AND ELECTRON 
BACKGROUNDS 

In order to directly compare our results to the observations 
made by HEAT, we must calculate the ratio of the local 
positron flux to the combined local positron and electron 
fluxes, called the "positron fraction" . In order to do this we 
require the background spectra for secondary positrons, pri- 
mary electrons and secondary electrons, where we propose 
that the primary positron component originates from the 
neutralino annihilations. In the following, we make use of the 
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param eterised fits to these spectra calculated bv lEdsio et alJ 
( 1998) and stated here (in units of GeV -1 cm -2 s _1 sr" 1 ), as 
follows: 



dE e + 
( d$ 



0.16e" 



\dE e 



1 + lle 0.9 +32e 2.15 

0.70e ' 7 



e ' / sec. e 



1 + llOe 1 - 5 + 600e 2 - 9 + 580e 4 - 2 



4.5e u 



dE e + ) scc e+ ~ 1 + 650e 2 - 3 + 1500e 4 - 2 



(4) 
(5) 
(6) 



again where e = E e + / (1 GeV). These equations agree with 
their respetive observational results to within 10-15% for the 
relevent energy intervals. (For a more detailed account of t he 
precision of these equations, refer to IStrong et alJ |2fi 
The positron fraction is then calculated using 
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5 DARK MATTER SUBSTRUCTURE 

The standard cosmological model assumes that all struc- 
ture in the universe originated from small amplitude quan- 
tum fluctuations during an epoch of inflationary expansion 
shortly after the big bang. The linear growth of the resulting 
density fluctuations is then completely determined by their 
initial power spectrum 



P(k) = k n , 



(8) 



where for n ^ 1, clumps are formed with a wide range of 
scales oc fc -1 . During the expansion of the universe, smaller 
clumps coalesce to form larger ones in a process of "bottom- 
up" hierarchical structure formation, in which the mass dis- 
tribution at any given redshift can be determined through 
the use of numerical simulations. 

In this paper we utilise the mass distribution of clumps 
deduced from the high-resolution simulations conducted by 
Diemand et al. (2005), of the form dnn/dlag(M/M ) oc 
(M/M )" 1 exp[-(A//Af cut _orf)~ 2/3 ]. Diemand et al. de- 
duced a lower mass cut-off M cu t_ fi ~ 8.03 x 1O~ 6 M0, 
however we shall also investigate the effect of using a much 
larger value, of order 10 6 Mq (similar to the typical mass res- 
olution of conventional large-scale numerical simulations), 
to simulate the destruction of lighter clumps by merging. 
We normalise the distribution so that the number density 
of clumps between 1O" 6 M and 1O _5 M is 500 pc~ 3 , with 
a halo-to-halo scatter factor of approximately 4, as stated 
by Diemand et al.. We also invoke an upper mass cut-off, 
M max ~ 1O 1O M , defined so that heavier clumps possess a 
number density of less than 1 per unit volume within the 
region in which positrons contribute to the local flux. 

The rate of dark matter annihilations within a clump 
crucially depends on its density profile. In this study we 



adopt the widely used Navarro, Frenk fc White (here after, 
NFW) density profile jNavarro. Frenk fc White! Jl99it ) 



p(r) 



Po 



c(r/r 200 )[l + c(r/r 20 o)] 2 ' 



(9) 



which is consistent with simulated results for the outer re- 
gions large-scale halos. The NFW profile also provides a 
good fit to those of the light est clumps produced in the 
simulations by Diemand et al. jDiemand et alJ J2005)). The 
maximum radius of a clump is defined here to be the radius, 
i~200, a t which its density is equal to 200 times the critical 
density. The concentration parameter c is a further degree 
of freedom, found to lie between 1.6 < c < 3 for the clumps 
in Diemand et al.'s simulations. Since the relationship bete- 
ween clump mass and concentration is complex, we approx- 
imate that all clumps have the same concentration and we 
compare results when using c = 1.6 and 3. 

Since the simulations conducted by Diemand et al. were 
terminated at z = 26, even though there is reason to believe 
that the shape of the mass distribution of clump s is similar 
today for large mass scales JDiemand et alJ J2005h ). we must 
also consider mass losses from tidal stripping by stars as 
clumps traverse the galactic halo. In this study we consider 
a scenario where the fractional mass loss experienced by 
each clump, over a given duration, is a mass-independent 
constant 1 — / and we compare results when using several 
different values of this parameter. 

The notion of a universal tidal stripping parameter is 
consistent with the calculation by Zhao et al., who demon- 
strates that 



111 4 oc-^ At, 

J Pcorc 



(10) 



where (p„) is the average stellar densiity along the clump or- 
bit, pcore is an effective "core" density of the clump and At 
is the time period considered. Sin ce the major i ty of clumps 
have similar orbital parameters (Zha o et alJ ||2p05) ) , and 
since we have adopted a universal clump profile, equation 
1101 1 implies that / is mass-independent. 



6 POSITRON SPECTRA FROM DARK 
MATTER SUBSTRUCTURE 

We consider four different neutralino models: 

First, a 50 GeV neutralino with an annihilation branch- 
ing ratio of 0.96 to bb and 0.04 to t + t~ . Such a particle could 
be gaugino-like or higgsino-like, since for masses below the 
gauge boson masses, these modes dominate for either case 
(designated as model 1). 

Second, we consider two cases for a 150 GeV neutralino. 
One which annihilates as described in model 1 (designated 
as model 2) , and another which annihillates entirely to gauge 
bosons, W + W~ or ZZ (designated as model 4). Such neu- 
tralinos are typically gaugino-like and higgsino-like respec- 
tively. 

Finally we consider heavy, 600 GeV neutralinos, which 
annihilate to bb with a ratio of 0.87 and to t + t~ or t + t~ the 
remaining time (designated as model 3). (We note that even 
though we do not explicitely calculate results for a model in- 
volving 600 GeV Higgsino-dominated neutralinos, which pri- 
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mx = 50 GeV 
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Figure 2. The calculated positron fraction as a function of 
positron energy (in GeV) , for a 50 GeV neutralino which anni- 
hilates to bb 96% of the time and 4% to t+t — (designated as 
model 1). The error bars displayed are for the 1994-95 (red) and 
2000 (blue) HEAT data. The solid line represents the spectra 
which best fits the data. The best-fit x 2 f° r this model is ap- 
proximately 5.4, (for 12 data points). Dashed lines represent the 
spectra corresponding to the 1 a fits to the observations, where 
the normalisation of the positron flux was considered to be a free 
paramter. The dot-dashed line represents the non-exotic back- 
ground contribution to the positron fraction. 



mx = 150 GeV (Gaugino) 
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Figure 3. The calculated positron fraction as a function of 
positron energy (in GeV), for a 150 GeV neutralino which an- 
nihilates to bb 96% of the time and 4% to t + t~ (designated as 
model 2). The best- fit x 2 for this model is 4.4, (for 12 data points). 
Otherwise, the same as in figure 2. 

marily decay to gauge bosons and partially to Higgs bosons, 
the resulting positron spectra per annihilation is likely to be 
very similar to that of model 3) . 

Although these models do not fully encompass the ex- 
tensive parameter space available to neutralinos at present, 
they do describe effective MSSM benchmarks. Furthermore, 
the relevant results for neutralinos with a mixture of the 
properties of those above can be inferred by interpolating 
between those presented. 

Figures F2l5l display the positron fraction as a function 
of positron energy, calculated for each of the four MSSM 
models outlined above. 

The thermally-averaged product of the low- velocity an- 
nihilation cross-section and relative speed of dark matter 
particles, (crv), was left as a free parameter and varied to 
fit our results to observations. In each figure, the solid line 
displays the positron fraction which fitted best to 1 a error 
bars of the 94-95 (red) and 2000 (blue) HEAT data. Dashed 
lines correspond to positron fractions where x 2 (for the 12 
data points), differs by unity from the corresponding best- 



mx = 600 GeV (Gaugino) 
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Figure 4. The calculated positron fraction as a function of 
positron energy (in GeV), for a 600 GeV neutralino which annihi- 
lates 87% to bb and 13% to t+t — or t+t - (designated as model 3). 
The best-fit x 2 f° r this model is 2.9, (for 12 data points). Other- 
wise, the same as in figure 2. 
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Figure 5. The calculated positron fraction as a function of 
positron energy (in GeV), for a 150 GeV neutralino which annihi- 
lates entirely to W + W~ or Z + Z~ (designated as model 4). The 
best-fit x 2 f° r this model is 3.3, (for 12 data points). Otherwise, 
the same as in figure 2. 

fit value (la results). Finally, the dot-dashed lines corre- 
sponds to the non-exotic background spectrum. The results 
are summarized in tables and [2] 



7 ASSESSMENT 

The results displayed in figures I2T51 indicate that neutralino 
annihilations within local dark matter substructure can 
give rise to positron fractions which correspond well with 
observations, with \ 2 between 2.9 and 5.4, for values of 
(av) ranging from 10 _29 cm 3 s _1 for 50 GeV neutralinos to 
2 x 10" 27 cm 3 s" 1 for 600 GeV neutralinos. 

Such values of (av) are at least an order of magnitude 
smaller than the canonical estimate of 3 x 10~ 26 cm 3 s -1 , 
determ ined from relic density calulations feertone et alJ 
(2005)). However, many conventional MSSM models exist in 
which this approximation is grossly violated and are easily 
able to accomodate the results obtained, where the num- 
ber of such models rapi dly increases for ligh ter neutralinos, 
especially near 50 GeV IIHooper et alJ ((2005)). 

As expected, we observe that the best-fit values of (av) 
increase as the dark matter clumps become increasingly 
stripped, and decrease as we increase their concentration, 
where the latter makes sense since a larger concentration 
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10-29 


Model 1 
1.0 
0.1 
0.01 


1.40 
2.30 
5.38 


0.55-2.05 
1.39-3.21 
3.41-7.36 


3.99 
6.61 
16.61 


2.51-5.48 
4.19-9.03 
10.54-22.68 


Model 2 
1.0 
0.1 
0.01 


2.99 
4.42 
10.08 


1.92-4.07 
2.88-5.95 
6.58-13.58 


7.46 
12.37 
31.09 


4.84-10.09 
8.10-16.65 
20.34-41.84 


Model 3 
1.0 
0.1 
0.01 


10.90 
18.04 
42.42 


6.47-15.34 
11.82-24.26 
28.99-55.85 


31.47 
52.05 
131.90 


21.36-41.58 
35.60-68.50 
90.21-173.59 


Model 4 
1.0 
0.1 
0.01 


3.81 
5.56 
12.52 


2.56-5.06 
3.76-7.36 
8.20-16.85 


8.73 
15.68 
39.82 


5.44-12.02 
10.44-20.92 
26.77-52.86 


Table 1. Values of the product (<?v), in cm 3 s 1 , which give rise 
to positron fractions which best-fit the HEAT observations, for 
the four MSSM models considered, for M cut _ off = 8.03 X 10~ 6 
and various values of c and /. One sigma error ranges are also 
displayed. 
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Model 1 
1.0 
0.1 
0.01 


6.52 
9.57 
21.90 


4.07-8.98 
6.05-13.08 
13.81-29.99 


16.14 
26.68 
67.69 


10.05-22.24 
16.39-36.97 
42.95-92.44 


Model 2 
1.0 
0.1 
0.01 


11.75 
17.13 
36.97 


7.63-15.87 
11.22-23.04 
22.43-51.52 


28.87 
48.47 
121.96 


18.48-39.26 
31.67-65.28 
79.81-164.10 


Model 3 
1.0 
0.1 
0.01 


47.73 
72.36 
163.78 


31.16-64.30 
48.95-95.74 
111.20-215.35 


122.08 
201.55 
508.81 


83.55-160.60 
137.95-265.14 
348.73-668.90 


Model 4 
1.0 
0.1 
0.01 


12.27 
20.85 
49.64 


6.94-17.61 
13.39-28.31 
33.57-65.71 


36.76 
60.95 
154.71 


24.69-48.83 
41.29-80.61 
104.77-204.65 



Table 2. Values of the product {crv), in cm 3 s — 1 , which give rise 
to positron fractions which best-fit the HEAT observations, for 
the four MSSM models considered, for M cut _ fj = 10 6 and vari- 
ous values of c and /. One sigma error ranges are also displayed. 
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implies that an increasing proportion of the clump mass is 
contained within a smaller radius. However, we observe that 
even with up to 99% stripping, and varying the clump con- 
centration over the entire range observed by Diemand et al., 
the best-fit values of {crv) do not change by more than an 
order of magnitude. 

We note that, of course, in reality the density of the 
clumps cannot diverge at their centres, as an NFW profile 
suggests, but are likely to possess unifor mly dense cores co n- 
taining of order 1% of their total mass (|Zhao et alJ J^005 ) ) . 
The stripping and eventual disintigration of these cores is 
a topic of much de bate and is currently being investigated 
iGreen et alJ (120061) '). We therefore do not display results for 
/ < 0.01, since it is unknown as to whether these cores can 
survive such severe levels of disruption. 

However, in table|2]we display results similar to those 
in tabled now for M cut _off = 10 6 Mq, corresponding to the 
typical mass resolution of simulations of the galactic halo, 
in order to simulate the destruction of lighter clumps. We 
observe that the best-fit values of (crv) increase by no more 
than an order of magnitude compared to their corresponding 
values in tableQ but are more consistent with the canonical 
value and consequently consistent with a larger number of 
MSSM models. 

We should also note that these results may be slightly 
modified when acknowledging flux contributions from tidal 
streams res u lting from material stripped from clumps 
iZhao et alJ J2005)). However, owing to the significantly 
lower densities of these streams, compared to the clumps 
from which they originate, it is very likely that such effects 
will not significantly alter our results. 



8 CONCLUSIONS 

For each of the four benchmark supersymmetric models con- 
sidered, we solved the diffusion-loss equation for the local 
positron flux resulting from neutralino annihilations within 
local dark matter substructures. We utilised the results from 
the unprecedentely high resolution simulation conducted by 
Diemand et al., who determined that the lightest clumps to 
form are of order an Earth mass and as large as 10 -2 pc. 
The mass distribution for clumps resulting from the simula- 
tion, which was terminated at z = 26, has been shown to be 
similar to be that expected today for large mass clumps. We 
assigned each clump an NFW density profile, with an as- 
sociated universal concentration parameter consistent with 
the numerical results. We investigated the effects on the flux 
when invoking different degrees of (universal) mass loss suf- 
fereed by each clump, in order to simulate the effects of 
tidal-stripping from stellar encounters. 

Our results indicate that neutralino annihilations 
within local dark matter substructure can give rise to 
positron fractions which correspond well with the observa- 
tions by HEAT, for values of {crv) ranging from 10 _29 cm 3 s _1 
for 50 GeV neutralinos to 2 x 10" 27 cm 3 s _1 for 600 GeV neu- 
tralinos. Despite that such {crv) are at least an order of mag- 
nitude smaller than the canonical value of 3 x 10 -26 cm 3 s _1 , 
determined from relic density calulations we indicate that 
many conventional MSSM models exist in which this ap- 
proximation is grossly violated and are easily able to ac- 
comodate our results obtained, where the number of such 
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models rapidly increases for lighter neutralinos, especially 
near 50 GeV. 

We observed that even with up to 99% stripping, and 
varying the clump concentration over the entire range ob- 
served by Diemand et al., the best-fit values of (ov) do not 
change by more than an order of magnitude. We also ob- 
served that even if we change the lower mass cut-off in our 
mass distribution from an Earth mass to 10 6 solar masses, 
the best-fit values of (av) increase by no more than an order 
of magnitude, but are consequently more consistent with a 
larger region of MSSM parameter space. 



APPENDIX: DERIVATION OF THE 
STEADY-STATE POSITRON FLUX 

Here we derive in detail the steady-state solution to the 
diffusion-loss equation given by 



d dn e 
didK 



K(E e+ ,r)V 



dE e+ 



b(E e+ ,r 



dn e + 



dE„ 



+ Q(E e+ ,r). (11) 



Our treatment is very similar to that described in 
lEdsio et al] <ll998T) . except that several unnecessary assump- 
tions made in that derivation are omitted here, making our 
solution slightly more general. 

We firstly re-arrange 1111 in terms of the dimensionless 
parameter u = 1/e = 1 GeV/iJ e + , and re- write the diffusion 
coefficient, given by as K — Koh(u), which gives 



1 d dn e+ V7 2dra e + r 2,, s 



Q(e(u),r), 
(12) 

where we have used equation J2J for the energy-loss rate, 
b. We now re-express 1121 in terms of the variable v, where 
h(u) = dv/du. Provided that h(u) is differentiable, we obtain 



d dn p 



= KoteV 



dv du du 
Now making the substitutions 



2 ^-te [u 2 h(u)] 1 Q(e(v),r). 



w(v) 



de 
dv 



= [u 2 h(u)] 1 , 



and 



F(v, r) 
equation i l-'it becomes 



dn e +(v,r) 
du(v) ' 



(13) 



(14) 



(15) 



■^-F(v,r)=KoT E V 2 F(v,r) 



■T E w(v)Q{e(v),r). (16) 



As stated in lEdsio et al] (Il99fif) . equation 1161 is analo- 
gous to an inhomogeneous heat equation. The spatial vari- 
ables are exactly analogous, while v is analogous to "time". 
Since v = J h(u)du — 3 a u + (1 — q)" 1 ^ 1 ^", we find that v is 
a monotonically increasing function for a < 1, and therefore 
is appropriately analogised to time. 

Considering the above, we can write the solution of the 



steady-state positron number density as an integral over the 
effective source term, — TEw{v)Q(e(v), r), of equation H16II . 
multiplied by a Green's function. 

We firstly solve for the free Green's function Gf re e, 
which satisfies 



Gfree (l 



KoteVG{ 
S(v — v')S(r - 



(17) 



which is identical to equation Ijlfip with an effective source 
term 5(v — v')8(r — r'), which is proportional to the source 
term for a monoenergetic point source, with v = v', located 
at r = r' . Taking the (spatial) Fourier transform of 1)17^ , we 
obtain 



Gfree (v ~ V , k) + K T E V 2 fc 2 Gfree (v - V , k) 

= (2ir)~ 3 ^ 2 5(v — u')exp(— ik ■ r'), 



(18) 

where Gfree is the Fourier transform of G{ lac . Multiplying 
each side of Ijl8fl by the integrating factor I = exp(KoTEk 2 v) 
and integrating over v, we obtain 



G{v-v',k) = {2TTy 3/2 oxp(-K TEk 2 (v-v')) 

xexp(— ik-r)0(v -v'), (19) 

where we have invoked the condition G{ lac = (and there- 
fore Gfree = 0) for v — (i.e. at infinite energy). Finally, we 
take the inverse Fourier transform of 1191 . to obtain 



Gfree(« - v, r - r') = (2iry 3/2 J G(v - v' , 
which evaluates to 



G bee (v-v',T-r') = [D(v,v')]~ 



exp 



k)exp(ik • r)d k, 
(20) 



[D(v,v')] 



(21) 



where D(v,v') = [4mKoTE{v — v')] 1 ^ 2 . 

As described in § 3, our chosen diffusion zone is an in- 
finite cylindrical slab of thickness 2L, where positrons lo- 
cated outside of this region simply free-stream. Therefore 
the steady-state positron distribution must therefore vanish 
at r = oo and \z\ > L. Ignoring contributions form sources 
outside of the diffusion zone, then from the form of the dif- 
fusion equation, the positron distribution will be zero for 
\z\ > L if it is zero at z = ±L. Therefore in order to obtain 
our desired positron distribution within the diffusion zone, 
we can use a Green's function which simply vanishes at the 
boundaries of our diffusi on zone . To do this we adopt the so- 
lution proposed bv lEdsio et al] (^998), which utilises a set 
of image charges, 



Vn=V 



2Ln + (~l)V 



(22) 
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to find the required Green's function, 



Considering the above, the source function is given by 



G 2L (v-v',r-r') = ]T (-l) n G bca (v - r - r'„). (23) 



n — — oo 



As described above, the steady-state positron distribution 
can then be expressed as the following integral 



— = te€~ / dv'wlv') 
Ot Jn 



x J d 3 r'G 2 L («(c) - v', r - r') Q(e(t/), r'). (24) 

Then using equation 1141 . we can conveniently re-express 
FM as 



9n 



= r E e 



de' 



x y d 3 r'G 2L («(c) - v'(e'), r - r') Q(e, r'). (25) 

Now we must construct an appropriate source func- 
tion, Q. To do this we only consider contributions from 
dark matter substructure, ignoring the contribution from 
the smooth component of the halo, which we calculate to be 
approximately 10 -4 times the substructure component. As 
described in § 4, we assume each clump to have an NFW pro- 
file, Pnfw(c, ?"20o(M)), with concentration, c, and maximum 
radius, r2oo(c, M), determined by c and the clump mass, M. 
The rate of production of positrons of energy e GeV, is then 



Q{e,r) 



(erf) d(j) 



/nfw (c,M) 



dn(r, M) 
dlogM 



dlogM, 



(29) 

where M m i n ~ 10 6 M Q and M max ~ 10 10 M Q , as discussed 
in § 4. 

Substituting expressions 1231 . 1291 and (12811 into equa- 
tion 1251 . evaluated at z — and r = Rq, we obtain the 
required (local) positron number density 



de 



(m x e) 2 



/Spw (c,M) 



dn D (M) 
dlogM 



dlogM 



Aim!,, 

de'{D[v(e),v'(e')]}- 3 ^ g 



Ti — — OC 



x dr r 



(-l) n z' + 2Ln 
D[v(e),v'(e')} 

0.86 



(r'/R)(l + (r'/R)) 2 



xexp — 



2rr 

x y dO'exp 



{D [v(e),v'(e')]} 2 
2rr'cos8' 



{D[v(e),v'(e')]y 



(30) 



m x de 



(aw) d<^ poofeM) ( CTW ) d 

' p (r )47rr dr = — / NFW (c,Af) 

m x de 

(26) 

where (cru) is the thermally-averaged, low-velocity annihi- 
lation cross-section multiplied by relative speed and d<f>/de 
is the number of positrons produced per annihilation per 
energy interval. 

We must now multiply equation 1261 by the correctly 
normalised density of clumps. As described in § 4 we utilise 
the clump distribution, dno{M) / dlogM , determined by Die- 
mand et al., normalised to a local clump density of 500 pc -3 , 
and like Diemand, we assume that this normalisation is pro- 
portional to the underlying halo density profile, phaio- For 
simplicity, we choose the halo profile to be a 'cylindrically' 
symmetric NFW profile 



Phalo(r) 



Po 



(r/R)(l + (r/R)) 2 



(27) 



where r is the cylindrical radial coordinate, R ~ 20kpc 
and po is a constant. The correctly normalised, spatially- 
dependent mass distribution of clumps is then 



dn(r,M) 
dlogM 



Phaio(r) dn D (M) 



Phaio (r = Re) dlogM 

0.86 dn D (M) 
(r/R)(l+(r/R)) 2 dlogM ' 



(28) 



where Rq ~ 8.5 kpc, is the approximate distance between 
the Earth and the galactic centre. 



The 0' integral evaluates to 27rJo, where 7o is a modified 
bessel function of the first kind and the z integral evaluates 
to the following sum of error functions, 



H(L, D) = 77 1/2 £>(erf (^) + ]T j erf 



+erf 



(4m + l)L 



(4m - 3)L 



D 



D 



2erf 



(4m - 1)L 



D 



(31) 



We are then left with integrals over e and r , 



' /nfw(c, M) // dlogM 



dn D (M) 



(™ x e) 2 J dlogM 

OO 

x J de' [D(v y v')]- 2 ^H(L,D(v y v')) 



x / dr' 



0.86 



o 

xexp 
x/ 



(l + (r'/R)) 2 

g| + r' 2 
[D(v,v')] 2 
2R Q r' 



,[D(v,v')] 2 ' ' 
which we evaluate numerically. 



(32) 
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Finally, using equation 13211 . we obtain the local 
positron flux resulting from neutralino annihilations to be 

V / local V / local 

where (5c is the speed of a positron of energy e. 
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